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Abstract 

We introduce and analyze a Nitsche-based domain decomposition method for the solu- 
tion of hypersingular integral equations. This method allows for discretizations with non- 
matching grids without the necessity of a Lagrangian multiplier, as opposed to the traditional 
mortar method. We prove its almost quasi-optimal convergence and underline the theory by 
a numerical experiment. 
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1 Introduction 

We propose and analyze the Nitsche method as a simple domain decomposition method for the 
solution of hypersingular boundary integral equations. In this context, simple means that (i) its 
implementation is not more difficult than a conforming approach and (ii) its numerical analysis 
avoids mathematical difficulties inherent to usual domain decomposition approaches. Still, a 
thorough analysis of our method, given in this paper, faces the problem of non-existence of a 
well-posed continuous counterpart for the discrete formulation. This is due to the low regularity 
of the underlying energy space. Main attraction of the Nitsche method, apart from its relative 
simplicity, is that it can maintain ellipticity and symmetry of the original problem. 

We study the hypersingular integral equation governing the Laplacian in IR^ exterior to 
an open surface, subject to a Neumann boundary condition. In principle, our domain decom- 
position approach is applicable to more realistic problems like linear elasticity and acoustics. 
Nevertheless, whereas a generalization to the Helmholtz equation is not difficult (it is a compact 
perturbation of the Laplace case) there are major difficulties in case of the operator governing 
the Lame equation. This remains an open problem. 
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For the solution of partial differential equations, domain decomposition is a classical strategy. 
It is used mainly for parallelization and the solution of linear systems. A variety of techniques 
exist, such as alternating Schwarz methods (see e.g. [26j). Of particular interest are methods 
that allow for non-matching meshes at the interface between sub-domains. They facilitate to a 
great extent mesh generation for complicated geometries. The so-called mortar method has been 
designed for this purpose [H [9] . It consists in introducing an unknown Lagrangian multiplier on 
the interface and adding interface conditions in a weak sense. For an analysis of the Laplacian 
in two and three space dimensions see [7J. This method transforms the original problem into a 
saddle-point structure, so that any numerical scheme requires a discrete inf-sup condition, i.e. 
compatibility between approximation spaces on sub-domains and the interface. 

An alternative to the mortar method is Nitsche's method, originally published in |25l [T]. and 
adapted in [6] to a domain decomposition framework. The interface condition is again treated 
weakly; not as an additional equation but like a penalization term in the (discrete) variational 
formulation. Other terms are added to the formulation to achieve consistency and ellipticity. 
Moreover, symmetry can be maintained for symmetric problems. As a result, Nitsche's method 
differs from classical penalization methods where consistency is lost [6]. 

In conclusion, main advantages of Nitsche's method are that 

1. no additional unknown is needed on the interface, 

2. no inf-sup condition must be satisfied among discrete spaces (except for the global ones, 
of course), and that 

3. discrete problems are elliptic and can be symmetric for symmetric problems, so that 

4. standard linear solvers can be used. 

Nitsche's method is closely related to the stabilized method of Barbosa & Hughes O S], 
which also circumvents the inf-sup compatibility condition that arises when a Dirichlet boundary 
condition is imposed weakly through Lagrangian multipliers. The connection between the two 
methods is established in |28j . 

In the context of partial differential equations, the Nitsche method has been applied suc- 
cessfully to a variety of problems such as linear elasticity [HJ [5], two- phase flows [27], and 
fluid-structure interaction [171 El H] • 

In the context of boundary integral equations and the use of non-matching grids or weakly 
imposed boundary (or interface) conditions, we only know of the results \15\ I18j. Both analyze 
a setting based on Lagrangian multipliers. The former reference provides the basic results 
like an integration-by-parts formula for the hypersingular integral operator, and analyzes the 
implementation of Dirichlet boundary conditions in a fractional order Sobolev space of order 
1/2. The latter reference proposes and analyzes the mortar domain decomposition approach for 
the hypersingular integral equation. An extreme case, the use of discontinuous basis functions 
for hypersingular operators, is studied in |20j . 

Let us also mention that there are several papers on domain decomposition involving bound- 
ary elements, e.g. [2lj where standard boundary elements are used for problems on sub-domains 
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of the PDE problem, and |301 ll9j which analyze domain decomposition for boundary elements in 
the construction of preconditioners. These papers do not deal with the problem of approximating 
functions (of fractional order Sobolev spaces) in a non-conforming way. 

In this paper we propose and analyze a Nitsche domain decomposition variant for the hy- 
persingular integral equation governing the Laplacian. Although this approach is simpler than 
mortar strategies in important aspects, as explained before, there are some non-trivial obstacles 
in its numerical analysis. Energy spaces of hypersingular operators are fractional order Sobolev 
spaces of order 1/2. These spaces form the natural basis for variational formulations. Now, do- 
main decomposition introduces interfaces where discontinuities arise. In the variational setting, 
these discontinuities are not well posed, simply because no well-defined trace operator exists. 
Therefore, we analyze the discrete Nitsche method without using a corresponding variational 
formulation. This is very much in the spirit of Strang's second lemma for non-conforming meth- 
ods. The difficulty of non-existence of a well-posed trace operator reappears in the analysis of 
the discrete problem. We deal with this problem by making use of a whole scale of Sobolev 
spaces (of higher regularity than 1/2) and by using inverse properties of discrete functions. The 
result is an almost quasi-optimal error estimate for the Nitsche method. Here, "almost" refers 
to perturbations which are only logarithmic in the mesh size. 

The rest of this paper is organized as follows. In J}2]we define some Sobolev spaces and our 
model problem. We also briefly recall the standard boundary element approximation. In ^ we 
introduce a domain decomposition (for simplicity only into two sub-domains; but this generalizes 
to more sub-domains in a straightforward way), the Nitsche-based discretization, and present 
our main result (Theorem 13. ip . Technical details and the proof of Theorem 13.11 are given in ^ 
In ^we present some numerical experiments that confirm the theoretical result. 

Throughout the article, we will use the symbols "<" and ">" in the usual sense. In short 
O'hiv) ^ bh{v) when there exists a constant C > independent of v, the mesh size h and a 
fractional Sobolev index e (if present), such that: ah{v) < Chh{v). 



2 Sobolev spaces and model problem 

First let us briefly define the needed Sobolev spaces. We consider standard Sobolev spaces where 
the following norms are used: For C IR" and < s < 1 we define 

II l|2 _ II ||2 , t |2 

W^Ww'iQ.) ■— IpIIl2(q) -I- 

with semi-norm 

f f f \u{x) -u(2/)P , , \i/2 

For a Lipschitz domain and < s < 1, the space i7*(il) is defined as the completion of C^{p) 
under the norm 

_ A |2 f ^ V/"^ 

{Mhs(U) + ]^ (dist(x,517))2^ 
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For s £ (0, 1/2), || • and || • ||H''(r2) are equivalent norms whereas for s G (1/2, 1) there holds 

i/*(r2) = Hq{Q,), the latter space being the completion of C^{Q) with norm in H^{Q). Also 
we note that functions from /7^(0) are continuously extendible by zero onto a larger domain. 
For all these results we refer to \22\ I16j. For s > the spaces H~^{il.) and H~^{il.) are the dual 
spaces of if*(r2) and ff*(0), respectively. 

Let r be a plane open surface with polygonal boundary. In the following we will identify P 
with a domain in IR^, thus referring to sub-domains of F rather than sub-surfaces. The boundary 
of F is denoted by dV. 

Our model problem is: For a given function f G L'^{T) find u G H^^^iV) such that 

wu{x) := [ <y)j^r^ '^^y = ^(^)' ^ ^ (^.i) 

47r dn^ Jy oriy \x - y\ 

Here, n is a normal unit vector on F, e.g. n = (0,0,1)"^. Note that W maps ^^/^(F) 
continuously onto i?~^/^(F) (see [12]). The variational formulation of (|2.ip is: Find u G H'^/'^iT) 
such that 

{Wu,v)T = U,v)r Mv^H^'^iV). (2.2) 

Here, (■, •)r denotes the duality pairing between H~^/^{T) and H^/'^{T). Throughout, this 
generic notation will be used for the L^-inner product as well as for other dualities, the domain 
being mentioned by the index. 

A standard boundary element method for the approximate solution of ()2.2p is to select a 
piecewise polynomial subspace Hh C H^^'^iT) and to define an approximant Uh G Hh by 

{Wuh,v)r = {f,v)r e Hh- (2.3) 

Such a scheme is known to converge quasi-optimally in the energy norm. In ^ we will com- 
pare such a conforming approximation with our proposed Nitsche approach and a Lagrangian 
multiplier variant. 

3 Discrete variational formulation with Nitsche coupling 

In this section, we introduce the Nitsche-based boundary element method for the approximate 
solution of problem (|2.2p . and present the main result. Theorem 13.11 

3.1 Some preliminaries 

We consider a decomposition of F into two non-intersecting polygonal sub-domains Fi and F2. 
The extension to an arbitrary number of sub-domains is straightforward. We will denote this 
partition of F as 

r:={Fi,F2}. 
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The interface between the sub-domains is denoted by 7 := f i n f 2. Throughout the paper, we 
will use the notation Vi for the restriction of a function v to a sub-domain T,. Also, as in [6], we 
will use the following notation for the jump on 7: 

[v]:= (vi-V2)\y. 

Corresponding to the decomposition of T, we will need product Sobolev spaces, e.g. 

with usual product norm. This notation (putting the decomposition T instead of F) is used 
generically, i.e. also for the spaces H^{T). We introduce the following inner product 

{v,w)r ■■= {vi,wi)r^ + {v2,W2)r2 

for v,w € L^iT){= L^iT)) and its extension by duality to H%T) x H-''{T). 

For 1/2 < s < 1, we introduce the following (semi-)norms, that are needed for the error 
analysis: 

( 2 

(3.1) 



w 



2 

lk"(r) -^X] l^l?^"(r») + IIMIIl2(7)' 



i=l 



where \v\fjs{Y-^^ is the Sobolev-Slobodeckij semi-norm as previously defined. The case s = 1/2 
will be used only for discrete functions where the jump across 7 is well defined. 

To introduce the discrete scheme, let us define regular, quasi-uniform meshes %, i = 1,2, of 
shape regular elements (quadrilaterals or triangles): Fj = UxeTi^- "^^^ maximum, respectively 
minimum, diameter of the elements of Ti is denoted by hi, respectively h^, and we define: 

/i := max {/ii, /12}, Zii := niin{/i^, /i2}. 

Throughout this paper we assume that h < 1. This is no restriction of generality and is just 
needed to simplify the writing of logarithmic terms. We introduce discrete spaces on sub-domains 
consisting of piecewise (bi)linear functions: 

Xh,i '■= {v G C^{Ti); v\k is a polynomial of degree one for all K ^Ti, flarnsri = 0}, 

for i = 1, 2. We define a global discrete space on F: 

Xh ■= Xh,i X Xh,2- 

Note that functions v G do satisfy the homogeneous boundary condition along dV but are in 
general discontinuous across the interface 7. Therefore Xh (t H^^'^{T), and this discrete space 
cannot be used directly for the discretization ()2.3p . Instead, we reformulate ()2.3p as a Nitsche 
variant so that Xh can be used to approximate the continuous problem (12. 2p . 
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3.2 Setting of the Nitsche-based domain decomposition 

For the setup of the Nitsche method let us introduce the fohowing surface differential operators: 

curlip := [dx2^,-dx^ip,0), curl (p := dx^^P2 - d^^ipi for 9? = (/72, V's)- 

The definitions of the surface curl operators are appropriate just for flat surfaces (as in our 
case) but can be extended to open and closed Lipschitz surfaces (see e.g. [I0l[15]). We define 
corresponding piecewise differential operators curlf and curlj- as follows: 

2 2 
curlr : = ^ (curlr; ifif, curlr (f:='^ (curlr; (fif, 

1=1 i=l 

where curlr; and curlr; refer to the restrictions of curl and curl, respectively, to Fj, and {■)^ 
indicates extension by zero to F. We made use of the notation introduced before ipi = (p\ri, 
(^j = (p\ri- Furthermore, we need the single layer potential operator V defined by: 

Vifix) [ ^ dS„ if € {H-y\T))\ X G F. 

4vr Jy \x-y\ 

We define the following bilinear form on x Xh: 

Ar{uh,Vh) ■= {V curlr Uh, curlr Vh)r 

+ ^{TiUh - T2Uh, [vh])^ + ^{[uh],TiVh - T2Vh)^ (3.2) 

where > and a € {—1,1} are numerical parameters. The operators Tj are defined as follows: 

TiV=[{V curlr v)\rrU]\-y 

for 2 = 1,2. Here, tj is the unit tangential vector on 5Fj (in mathematically positive orientation 
when identifying Fj with a subset of IR^ which is compatible with the identification of F as a 
subset of IR^). Note that TiV is not well defined for v G H^/'^(T) in general since there is no 
well-defined trace from H^^^{Ti) to 5Fj. 

The Nitsche-based boundary element method associated to problem ()2.2p then reads: 

Find Uh € X^ such that 

Ar{uh,vh) = {f,vh)r (3.3) 

for all Vh G X^- 

Remark 3.1. For any function u £ H^{T) (s > 1/2), in particular for the solution of (12.10 . 
there holds {[u],[v]).y = {[u],Tiv — T2v)^ = for sufficiently smooth v. Therefore the terms 
([u/i], [v/i])^ and {[uh],TiVh — T2Vh)'y are not required for consistency of (13. 3p . However, the 
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additional term ^{[ufi],TiVfi — T2Vfi)'y in the Nitsche-based formulation is of interest for two 
reasons |3 Remark 2.11]. First, for a = 1, the bilinear form A'j-{-,-) becomes symmetric, as in 
the standard case (j2.2p . This allows in particular to make use of fast linear solvers for symmetric 
matrices. Also, for a = —1, symmetry is lost, but we recover ellipticity of A-i-{-, •) for any value 
of the parameter u > (see Lemma (i))- In- fo-ct, any value of a (including a = 0) can be 
chosen though only values —1 and 1 lead to interesting particular cases. 

The main result of this paper is: 

Theorem 3.1. Let u € H^'{T) with r S (1/2,1) be the solution of ()2.ip . In the case a = —1 
let u > be arbitrary and in the case a = 1 let v > Ci \ log /ip for a sufficiently large constant 
Ci > 0. Then, the discrete problem (|3.3|) is uniquely solvable and there exists a constant C2 > 0, 
depending on v and r, but not on u and the actual mesh, such that there holds the error estimate 

\\U - Uh\\ ^1/2 < C2|l0g/l|^/^ft,''"^/^||'u|lj:/r(r). 

A proof of this result will be given at the end of Section [H 

Remark 3.2. It is known that u € H^{r) for any r < 1, see, e.g., [29]. Using this regularity. 
Theorem 13.11 proves a convergence which is close to 0(/i^/^), the optimal one. The reduction 
to h"^^^^^ for any r < 1 is due to the assumed regularity in standard Sobolev spaces, and not a 
sub-optimality of the method. On the other hand, the logarithmic perturbation \ log/i^/^ is due to 
the Nitsche coupling, and is also present in non- conforming approaches (the same exponent 3/2 
appears in the Lagrangian multiplier approach [15], and in the mortar coupling \T8\ the exponent 
is 2). It is unknown whether these logarithmic terms in the upper bounds are optimal. 

4 Technical results and the proof of the main theorem 

In N4.ll we present some preliminary results and lemmas. In N4.2I we then prove the consistency 
of the method, using an integration-by-parts formula coming from [15^ I18j. Discrete continuity 
and discrete ellipticity are studied in We conclude with the proof of the main theorem in 

The steps followed in the error analysis are quite similar to those of the analysis of a Nitsche- 
based method for finite elements (see e.g. |6]). The main difficulty in the case of boundary 
elements consist in the non-existence of a well-posed variational Nitsche formulation. Error 
estimates are wanted in spaces related to H^^'^{T) where no well-defined trace operator exists. 
Therefore, the numerical analysis of ()3.3p makes use of a whole family of Sobolev spaces 
with r close to 1/2. Additional difficulty in our case is that operators are non-local in contrast 
to the finite element setting. In opposition to the mortar boundary element method [18], no 
inf-sup condition needs to be checked since no Lagrangian multipliers are introduced. 
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4.1 Preliminary results 

We first introduce the following spaces for the definition of the single layer potential operator 
V (see [H): 

hT\t) :={^G(F^-i(r))^ Vn = 0}, ^^^^ 
Hl{T) :={^G(if^(r))^ ^■n = 0}, 

where < s < 1 and the normal vector n has been defined previously. We will make use of the 
continuity (see [12]): 

V : hI~\t) ^ Hl{T), 0<s<l. (4.2) 
Lemma 4.1. For i = 1,2 there holds 

\\Tiv\\L2(^)< (s-l/2)-V2 ||curlr^;||^.-i(r) yvGH^T), 1/2 < s < 1, (4.3) 

\\Tiv\\m^)< (s-l/2)-3/2 1^1^^^^^ l/2<s<l, (4.4) 
\\TiVh\\L^{^) < [log/ip/^ \vh\Hi/2(j-) Mvh^Xh- (4.5) 

Proof. Let v G H^{T), with 1/2 < s < 1. We use the trace theorem |15|. Lemma 4.3] and the 
continuity of V (j4.2p to bound 



This proves (j4.3p . Further, by using the equivalence of the H^~^(Ti) and (Fj) norms for 

|s - 1| < 1/2 [19i Lemma 5] and the boundedness of curlr; : i^'*(Fj) i3"^~^(Fi) [181 Lemma 
3.4]), we obtain 



2 

|2 



(r»)- 



4 = 1 

1 ^ 1 ^ 

^ (7317^ E II -"rir. -llLr^(ro ^ (73T7^ E 11-11?^^ 

Combining these two estimates and using a quotient space argument proves ()4.4p . 
Now picking S X/^, we have: 

WTiVhU^-y) < (s - l/2)"^/^|7;/j|j^s(7-) 

using (14. 4p and then the inverse property \v\H^{r,) ^ Zi^^^^^l^'l//i/2(r,) (see, e.g., [191 Lemma 4] 
together with a quotient space argument). With the choice s = 1/2 + |log/i|~^, this proves 
[3D. □ 
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4.2 Consistency of the Nitsche formulation 

In this part, we show that the Nitsche formulation (j3.3p for the hypersingular operator is consis- 
tent, a classical result for the Nitsche method in the standard case (see e.g. [6]). One difficulty 
here is that the boundary operator Tj is not well defined for v € H^/'^{T). Nevertheless, we can 
take advantage of previous results proven in p5] . 

First, we need to start from an appropriate integration- by-parts formula for the hypersingular 
operator. For the convenience of the reader we recall the setting from [15^ 118]. For a smooth 
scalar function v and a smooth tangential vector field (p, integration by parts on Fj gives 

{(f ■ ti,Vi)er^ = (curlr; ip,Vi)r, - (curlr; Vi,(p)ri, 

for i = 1,2. We apply this formula to (p = curlti)[ri, so that: 

((ycurln)lr, ■ti,Vi)Qr^ = (curlri(Fcurlu),?jj)r, - (curlr; Uj, 1^ cur In) r^. 

Recalling the definition of Tj, and using a function v that vanishes on dV, we obtain 

{TiU,Vi)y = (curlr; curl u),7;j)r, - (curlr; -Uj, V curl n)ri. (4.6) 

Let us recall the following lemma from |18) (Lemma 3.5): 

Lemma 4.2. For u G H^/'^{T) with Wu = / G -^^(r); the equation (|4.6p defines TiU = 
(y curlu)[r, • t, G H-^i-f), with < s < 1/2. 

As a result, we can state: 

Lemma 4.3. Let > and |(t[ = 1. Then, the Nitsche formulation is consistent, i.e. the 
solution u of (12. Ij) (Wu = / G L^(F)j solves the discrete setting (13. Sh . 

^t{u, Vh) = if, VhW Vwft G Xh. 

Proof. Let u G #^/^(F) be the solution of (f2?T]) . It is weh known that u G F'*(F) for any s < 1, 
see, e.g., [29], so that the trace of n on 7 is well defined and [u] = 0. 

By Lemma [421 Tiu and T2U G H'^^^j) (0 < s < 1/2). Moreover, since [(ycurlu)jri ~ 
(y curlu)Ir2]|7 = and ti = — 12 on 7, there holds 

Tiu + T2U = on 7. (4.7) 

We obtain for Vh G 

ATiu,Vh) ={V curlr u, curlr Vh)T ~ ^2^, [vh])-y + ^{[u],TiVh - T2Vh)^ 

=(y curlr li, curlr t'/i)r +\{Tiu - T2U, [fh])^. 
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Using ()4.7p to write TiU = \TiU — \TjU {i ^ j), and rearranging terms, we obtain 



-^{Tiu - T2U, [vh])^ = '^{TiU,Vh,i) 



7 



i=l 



SO that, together with the integration- by-parts formula (j4.6p 



2 2 



Ariu,vh) = [{V curlr u, curlr; Vh,i)r, + {TiU,Vh,i)'y] = (curlr; curlr tt),^fe,i)ri 




=(curlr(y curlr n),f;/,)r = {Wu,Vh)r = {f,Vh)r, 



Here, we have also used the relation 



Wu = curlr V curlr u, 



see \23\ [23] . This proves the lemma. 



□ 



4.3 Discrete ellipticity and continuity 

Main advantage of the Nitsche method is that it yields elliptic bilinear forms in the case of 
elliptic problems. In the boundary element setting, we do not have an appropriate variational 
formulation. Nevertheless, discrete ellipticity is still achievable. This is contents of the first 
lemma. Afterwards, we briefly state discrete continuity without giving a bound for the continuity 
constant. This bound is studied in more detail in the proof of the main theorem in ^4.4[ 



Lemma 4.4. (i) Let a = —1. For all u > there exists a constant C{v) > such that 

Ar{vh,Vh) > C{v) \\vh\\\i/2,^. G Xh 



(ii) Let a = 1. There exists a constant Ci > such that, if u > Ci |log/ip, then there exists a 
constant C2 > independent of v such that 



and 




Ar{vh,Vh) > C2 \\vh\\'ti/2 



Vvh G Xh 



and 




Proof, (i) Case a = —1. 
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Let Vh € Xh- By the ellipticity of V and [151 Lemma 4.1] there holds 

2 

{Vcur\rVh,curlrVh)T > || curlr Vh|[^-i/2^j.^ > J]] || curlr; Wh||^_i/2^j.^^ 

2 

~ X] l'"'»lj/i/2(r,) = I ^^/i 1^1/2 (T)- 

i=l 

This proves that 

Ar{vh,Vh) = {VcurlrVh,curlrVh)T + '^{[vh],['^h])'y 
> \vh\Hi/2(r) + IIMIIl2(7) = 11'^'^ 1 1 Hi/2 (T) 
(which is the first assertion) and also 

Ar{vh,vh) > II curlr^h 11^-1/2^^^ + ||K]|li2(^). 

Both estimates together prove the second assertion, 
(ii) In the case o" = 1 we obtain for Vh € Xh 

Arivh,Vh) = {V curlrVh,curlrVh)T + J^iivh], [vkD-r + (^i^/i ~ '^'^Vh, [vh])'y 
Using the Cauchy-Schwarz and Young's inequahties, and (j4.4p . we bound 

{TiVh - T2Vh, K])7 ^ S\\TiVh - T2Vh\\l2(^^) + ^||N||i2(^) 

{s -472)3 ^ > 0- 



< 



Combination of these two relations with (j4.8p . and making use of the inverse property |w|ffs(r) ^ 
/ji/2-s|^|^^^^^^^ (see, e.g., [19< Lemma 4] together with a quotient space argument), yields 

AT{vh,Vh) > (l - SC, (,~\/2)3 ) i^^lHV2(r) + ('^ - ^)IINIli2(^) V5>0 
for two unknown constants Ci,C2 > 0. Selecting 

s = 1(1 + |log/i|-^) and S = ^\logh\-^ for c G (0, 1) 

this yields 

Ar{Vh,Vh) > bft|^l/2(7-) + ||Mlli2(^) = 1 1 1 1 5^1/2 

for v > ^ + c = ^^^^1 log /ip + c. This proves the first estimate in (ii). As in the case a = — 1, 
and using (j4.3p in addition to (|4.4p . one proves the second assertion under the same condition 
on u. □ 
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Lemma 4.5. Let > and \(t\ = 1. The bilinear form Aj- is continuous: 

Ar{vh,Wh) < C{iy,h)\\vh\\ ^1/2 ^^^\\wh\\ ^1/2^^^ '^Vh,Wh G Xh, 
with C{i',h) > a number that depends on v and on the mesh parameter h. 

Proof. This estimate follows by using the mapping properties of the involved operators V , curlT-, 
Ti and inverse properties of discrete functions. 

□ 

4.4 Proof of the main theorem 

By Lemma [4.4l and under the stated assumptions, the bilinear form Aj- is elliptic. Moreover, by 
Lemma 14.51 this bilinear form is also continuous on (with bound depending on the mesh) so 
that problem (|3.3|) has a unique solution. It remains to bound the error. To this end we follow 
the lines of a Strang estimate for non-conforming methods. By Lemma 14.41 there holds for any 

Vh G Xh 

\\U — Uh\\^l/2^^ < \\U — Vh\\j^l/2^^,^ + \\Uh — Vh\\j^l/2^^^ 

^11 II , ^riuh - vh,wh) 

< n-Wfe 1/2 + sup -| 1^ —- — . 4.9 

^* w,(,x^\{0} W curlr WhW ^- V2(p) + II [wh] ||l2(7) 

Now, by the consistency (see Lemma we obtain A'j-{uh — Vh,Wh) = A-j-{u — Vh,Wh) so that 
we continue bounding (using duality estimates and the continuity of V (j4.2p ) 

Ariuh - Vh,Wh) = (y curlr(u - Vh), curlr Wh)T + ~ '"h) - T2{u - Vh), [wh])^ 

+ ^{[U - Vh],TiWh - T2Wh)-y + 1^{[U - Vh], [Wh])-y 

< II curlr (-u -Vh)\\ ^-i/2(p) II curlr WhW ^-Va^p^ 

2 

+ ^[\\Tiiu - Wfe)||L2(7)IIK]||L2(7) + \\TiWh\\L^'y)\\[u - Vh]\\L2(j)^ 
i=l ' 

+ ll[^ - ^'/i]IIl2(7)||[u'/i]||l2(7). (4.10) 

We bound the terms on the right-hand side. 

In the following let s be a small positive number. Using \y%\ Lemma 5] and the continuity 
of curlp; : H^^^^'^^Vi) — > ll\ ""^'^^(Lj), together with a quotient space argument, yields 

2 

II curlr(u - w,,)f /2 < V II curlr(n - Vh)f--^i^ 

tit (i) -"t 

1 = 1 

2 , 2 



< ^ I] II curlr(^x - ^'/.)ll^pi/2(p^) ~ ^2 1^ - ^^^l^«+i/2(r,)' 

i=\ i=l 



(4.11) 
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\TiWh\\l2(^-) < -\\ curlrWhf^,_^/2 < — 2ji| curlr Whf^_y2,^,- (4.13) 



and estimate (|4.4p proves 

\\Ti{u - Vh)\\L2(^y) < - Vh\Hs+l/2(^r)- (4-12) 

Eventually, by (|4.3p and the inverse property (see [I9l Lemma 4]), 

1|| T l|2 < 1 

- curlr^fft -s-i/2 — n 

Combination of (j4.10p - (j4.13p proves that for any Wh G \ {0} and any small s > there holds 

Arjuh - Vh,Wh) _3/2/| I , 111 \ 

II curlru;^ 11^-1/2^^^ + II [u'ft]|Ii2(^) V y' > "-"J 

By a standard approximation result we have that, for r € (1/2; 1), 

Jnf_^(k - ''^/ilH''+i/2(r) + - Vh]\\L'^{^)) < /i''"'"^/^||'u|lj/r(r), 
so that referring to (I4.9P this proves that 

Ik - ^i/x||^i/2(^^ < h-'^'s-^l'^h'-^/^\\u\\Hr^^)■ 
Selecting s = | \ogh\^^ this proves Theorem 13. 1[ 

5 Numerical results 

We consider the model problem ()2.ip with F = (0, 1) x (0, 1) and / = 1. For the sake of simplicity, 
we only deal with the case of one sub-domain and in where the homogeneous Dirichlet condition 
on the boundary 7 = dV (implicitly present in the energy space H^f^iT)) is imposed weakly (in 
the discrete case) through a Nitsche formulation. 

This situation is identical to the one described in [151 Section V] where a Lagrangian mul- 
tiplier is used to impose the homogeneous boundary condition. Below we compare numerical 
results from both methods (Figure [5.5p . 

To obtain the Nitsche formulation of this problem, we formally extend n by onto IR^ and 
decompose IR^ into Fi = F and F2 = IR^ \F. The extension of n by is continuous in i/^/^(IR^) 
since u € H'^/'^iT). As a result, the Nitsche formulation is a particular case of the one studied 
in this paper, and the corresponding bilinear form is obtained from p.2p by using that jumps 
across 7 = dV are identical to traces on 7 (taking the exact approximation of the solution 
exterior to F). 

We use uniform meshes 7h on F which consist of squares of side-length h. The discrete 
spaces Xh are made of continuous piecewise bilinear polynomials on Th- Then the Nitsche-based 
formulation reads: Find Uh G Xh such that 

(ycurlr'Uh,curlr^;h)r + {Tuh,Vh)-y + cr{uh,Tvh)j + v{uh,Vh)'y = {f,Vh)r (5.1) 
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for all Vh G Xh. Here, the operator T is defined by 

Tv := t • y curlr 

with t being the tangential unit vector along 7. Since the exact solution u of ()2.ip is unknown, 
the error 

II l|2 _ I |2 I II ||2 

cannot be computed directly (note that u = on 7). Instead, we approximate an upper bound 
to the semi-norm \u — n/i|j|^i/2(-r) as follows. 
First, note that there holds 

\u - Uh\1ji/2(Y) < curlr(n - u/,),curlr(n - Uh))r, 



due to ellipticity of V. Taking into account that u is solution of (j2.ip and Uh is solution of (jS.ip . 
we find: 

\u - ^i/i|^i/2(r) < curlr n, curlr u)r + (T^ curlr Uh, curlr n/j)r - 2(y curlr u, curlr n/i)r 

= {Wu, u)r + (/, n/j)r - {Tuh, Uh)^ - a{uh,Tuh)-y 
- u{uh, Uh)'y - 2{Wu, Uh)r + ^Tu, u^)^ 

}y — 1^ \^hi "-hl^ T- U-, "'ft/7- 



< {Wu, u)t - if, Uh)v - (1 + (y){Tuh, Uh)j - v{uh,Uh)^ + 2{Tu, Uh)-y. 



Then, from {uh,Uh)'y > 0, the Cauchy-Schwarz inequality and Lemma HT] (inequality (I4.5p ). we 
obtain 

\u - -u/il^i/a^p) < {Wu,u)r - U,Uh)v + \\Tuh\\ 12 (^i^\\uh\\ 121^^1) + ||rti||^2(^)||-u,i||i2(^) 

< {Wu,u)t - {f,uh)r + (^|log/i|^/^j'U/i[^i/2(r) + \\Tu\\l2^^)^ ||u/,||i2(^). 

Note that for this specific problem, Tu G L'^il)- Furthermore, since the method is stable, 
|u/i|j:^i/2(-p) is bounded independently of h. This proves that 

\u-Uh\Hi/2(Y) < \{Wu,u)t - {f,Uh)v\^^'^ + I log /l I ^/"^ 1 1 U/, 1 1 ^''2^^^^. 
11I/2 

Moreover, since |p/i||l2(7) ^ ||^ih|l2,2(^)) there also holds 

h-^^/l|lj^l/2(p) < \{Wu,u)t - {f,Uh)T\^''^ +\\0^hf''^\\Uh\\^l2(^^y 

1 /2 

The terms {f,Uh)r and ||tt/i||L2(^) are easy to compute. The energy norm {Wu,u)j. of u can 
be approximated through extrapolation, denoted by || tijlea; in the following, see |13j . Therefore, 

\\H\l-{f,nh)r\'^'+ Ilog/ip/4h,|I^/4))/|ln||ex 
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is a computable and reasonable measure for an upper bound of the error \\u — Uh\\^i/2^^s^ nor- 
malized by ||^||j:^i/2(-r)- Below we present numerical results for the two contributions 

lll^llex - {f,Uh)r\^^^/\\u\\e^ (5.2) 
(referred to as "//i/^" error in the figures) and 

ex 

(5.3) 

(referred to as "L^" error). 

We first consider some tests in the skew-symmetric case {a = —1), for different values of z^. 
The corresponding results are given in Figures 15.11 and 15.21 A double logarithmic scale is chosen 
and the errors (j5.2p are plotted versus the dimension of the discrete space X^. Figure [5TT] presents 
results for the term (15. 2p of the error and indicates that the Nitsche-based method converges for 
all the tested values of z^, with a logarithmic perturbation of the convergence, as expected by 
the theory (Theorem 13. ip . As a consequence, the convergence is asymptotically a bit slower (by 
a factor of | log/i|) than in the case of the conforming BEM. The latter method converges like 
0(/i^/^), and for comparison we have given the curve | loghlh^^"^ as well (with a constant factor 
for adjustment). For all studied values of u the curves exhibit the same asymptotic convergence 
order, though their initial behavior differ. In particular, for v >2, a minimum is reached quickly, 
after which the asymptotic behavior is recovered. Apparently, for any particular mesh, there 
simply is an optimal value of z^ for term (15. 2p . 

Figure [5^2] shows that the other part of the error (given by ()5.3p ) also behaves as predicted. 
All the curves are parallel and the parameter does not seem to have a great influence, except 
for shifting the curves which corresponds to multiplication of the error by a constant. 

Next we study the symmetric case {a = 1). The corresponding results are given in Figures[5]3] 
and 15.41 As expected, if the value of v is not sufficiently large, the method does not converge 
(see the curve for = 1 in Figures 15.31 and l5.4p . Indeed, if v is too small, discrete ellipticity of 
Aj-{-, •) cannot be guaranteed (see Lemma l4.4p . Taking higher values of ly ensures convergence 
of the method. In particular, if is not taken as a constant but a power of | log h\ the asymptotic 
behavior improves. For u = \ log/ip, the behavior of the conforming BEM method is recovered, 
with quasi-optimal convergence. Note that theoretically, a sufficient condition to guarantee 
discrete ellipticity and convergence (in the symmetric case) is v > \logh\^ (cf. Lemma 14.41 and 
Theorem 13. ip . The same conclusions hold when one looks at the L^(7)-error (j5.3p in Figure [531 

In Figure 15.51 we compare the Nitsche method (symmetric and skew-symmetric versions) 
with the Lagrangian multiplier-based method [18]. One observes that the symmetric Nitsche 
method and the Lagrangian multiplier method have the same asymptotic convergence, which is 
quasi-optimal (without logarithmic perturbation) in this example. The skew-symmetric method, 
on the other hand, remains almost quasi-optimal, i.e. with logarithmic perturbation. 

Concluding, the numerical experiments are in good agreement with the theory, and illustrate 
the applicability of the Nitsche-based domain decomposition method for hypersingular integral 
equations, e.g. as a possible alternative to a Lagrangian multiplier approach which requires 
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|log h|*h'^{1/2} - \ F' 

conforming BEM — i — ^ ? 

Nitsche H"{1/2}, nu=1 — -<— - ■ / 

Nitsche H''{1/2}, nu=2 — ■ \ ^ 

Nitsche H»{1/2), nu=5 — »— / 

Nitsche H"{1/2}, nu=10 — a— \ I 

J , , , ii , I 

10 100 1000 10000 

dim(X_h) 



Figure 5.1: Skew-symmetric Nitsche method {a = —1): relative error curves (upper bound 
(|5.2|) ). Comparison with conforming BEM. 



an additional unknown and destroys ellipticity. In particular, the symmetric case seems to be 
more appealing due to its competitive convergence for large values of u, and since it maintains 
symmetry. 
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Figure 5.2: Skew-symmetric Nitsche method (cr 
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Figure 5.3: Symmetric Nitsche method (cr = 1): relative error curves (upper bound (|5.2 
Comparison with conforming BEM. 
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Figure 5.4: Symmetric Nitsche method (cr 
except for = 1, and (j5.3p ). 
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Figure 5.5: Comparing Nitsche and Lagrangian multipher methods: relative error curves (upper 
bounds ()5.2p and ()5.3p for all cases). 
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